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ABSTRACT 

Context. Observations of the Sunyaev-Zerdovich effect (SZE) from galaxy clusters are emerging as a powerful tool in cosmology. 
Besides large cluster surveys, resolved SZE images of individual clusters can shed light on the physics of the intra-cluster medium 
(ICM) and allow accurate measurements of the cluster gas and total masses. 

Aims. We used the APEX-SZ and LABOCA bolometer cameras on the APEX telescope to map both the decrement of the SZE at 
150 GHz and the increment at 345 GHz toward the rich and X-ray luminous galaxy cluster Abell 2163 at redshift 0.203. The SZE 
images were used, in conjunction with archival XMM-Newton X-ray data, to model the radial density and temperature distribution of 
the ICM, as well as to derive the gas mass fraction in the cluster under the assumption of hydrostatic equilibrium. 
Methods. We describe the data analysis techniques developed to extract the faint and extended SZE signal. We used the isothermal f} 
model to fit the SZE decrement/increment radial profiles. We performed a simple, non-parametric de-projection of the radial density 
and temperature profiles, in conjunction with X-ray data, under the simplifying assumption of spherical symmetry. We combined 
the peak SZE signals derived in this paper with published SZE measurements of this cluster to derive the cluster line-of-sight bulk 
velocity and the central Comptonization, using priors on the ICM temperature. 

Results. We find that the best-fit isothermal model to the SZE data is consistent with the ICM properties implied by the X-ray data, 
particularly inside the central 1 Mpc radius. Inside a radius of ~1500 kpc from the cluster center, the mean gas temperature derived 
from our SZE/X-ray joint analysis is 10.4 ±1.4 keV. The error budget for the derived temperature profile is dominated by statistical 
errors in the 150 GHz SZE image. From the isothermal analysis combined with previously published data, we find a line-of-sight 
peculiar velocity consistent with zero; v,. = -140 ± 460 km/s, and a central Comptonization yo = 3.42 ± 0.32 x 10"* for Abell 2163. 
Conclusions. Although the assumptions of hydrostatic equilibrium and spherical symmetry may not be optimal for this complex 
system, the results obtained under these assumptions are consistent with X-ray and weak-lensing measurements. This shows the 
applicability of the simple joint SZE and X-ray de-projection technique described in this paper for clusters with a wide range of 
dynamical states. 

Key words. Galaxies: clusters: individual: Abell 2163 - Cosmic microwave background - Cosmology: observations 



1. Introduction 

The Sunyaev-Zel'dovich effect (SZE, Sunyaev & Zel'dovich 
1970, Birkinshaw 1999) provides a powerful probe of the large- 
scale structure in the Universe by imprinting the thermal energy 
of clusters of galaxies on the cosmic microwave background 
(CMB). The inverse Compton scattering of CMB photons by the 
hot intra-cluster medium (ICM) has a characteristic frequency 
dependence and a power spectrum that is very different from 
that of the primary CMB anisotropies. 



* e-mail: mnord@astro . uni -bonn . de 



A number of blind SZE imaging surveys are in progress (e.g. 
Ruhl et al. 2004, Kosowsky 2006) with the aim of tracing the 
large-scale structure of the Universe through the detection of 
galaxy clusters, and the first blind SZE detections have recently 
been reported (Staniszewski et al. 2009). 

With sufficiently high resolution, the SZE can be used to map 
the pressure structure of individual clusters. Once combined with 
information on the X-ray surface brightness, the SZE thus pro- 
vides constraints on gas temperatures and total mass distribu- 
tions inside clusters. Such constraints are free of the potential 
biases of X-ray spectroscopy. A further comparison with weak- 
lensing maps yields insight into the dynamical state of the cluster 
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by checking the validity of the hydrostatic equilibrium condition. 
Joint de-projection methods using X-ray and SZE maps have 
been proposed (e.g. Lee & Suto 2004, Puchwein & Bartelmann 
2006), but so far their implementation using real SZE data has 
been limited to unresolved or barely resolved SZE maps (e.g. 
De Filippis et al. 2005). Kitayama et al. (2004) presented the 
first tempertaure deprojection of a cluster using a combined SZE 
and X-ray analysis. However, this analysis was limited to using 
a parametric model for the SZE signal. One main objective of 
this paper is to show the potential of such a de-projection anal- 
ysis using high signal-to-noise resolved SZE maps without the 
limitations of parametric models. 

Abell 2163 is a hot. X-ray luminous galaxy cluster at z = 
0.203, with a mean X-ray temperature in the central region of 
Jx = 12![ J keV (Markevitch & Vikhlinin 2001). It has a rel- 
atively large angular extent, with a virial radius estimated from 
weak-lensing mass modeling to be on the order of 15 arcminutes 
(Radovich et al. 2008), making the primary CMB anisotropics 
a major source of systematic uncertainty in the SZE measure- 
ments. Detailed X-ray observations suggest that Abell 2163 is 
a merger system (Elbaz et al. 1995, Govoni et al. 2004), with 
asymmetric X-ray temperature structures and strong radio ha- 
los. The merger scenario is also supported by optical observa- 
tions, most notably from the presence of two bright cD galax- 
ies (BCGs; Maurogordato et al. 2008). The cluster has been ob- 
served in the SZE at 30 GHz with OVRO/BIMA (Reese et al. 
2002), and at 142, 217, and 268 GHz using SuZIE (Holzapfel et 
al. 1997a,b). Using relativistic corrections to the SZE, Hansen et 
al. (2002) used these data to directly constrain the mean ICM gas 
temperature in Abell 2163, albeit with very large uncertainties. 

During the commissioning observation of the APEX-SZ 
bolometer camera (Schwan et al. 2003, Dobbs et al. 2006) in 
2007, we detected Abell 2163 with high significance (12cr) at 
150 GHz with arcminute resolution. We observed the cluster in 
the SZE increment at 345 GHz using the LABOCA bolometer 
camera (Kreysa et al. 2003, Siringo et al. 2009) on APEX, pro- 
viding images with an angular resolution of 19.5". 

With a map larger than 12' across, the present LABOCA data 
on Abell 2163 provide the first large-area imaging of a galaxy 
cluster at sub-mm wavelengths. High-resolution SZE increment 
imaging is important in view of planned future observations of 
clusters at sub-mm wavelengths, and for controlling potential 
foregrounds. However, due to limitations in constraining the SZ 
emission at large radii, the use of the sub-mm data in this paper 
is limited to the application of the multi-frequency SZE mea- 
surement to model the SZE spectrum and derive constraints on 
the thermal and kinematic SZE components. 

In sections |2] and [3] we describe the observations and data 
reduction algorithms used to obtain maps at 150 and 345 GHz. 
We also describe the analysis of complementary XMM-Newton 
data. In section |4] a simple isothermal model of the cluster 
gas density distribution is presented. A joint SZE/X-ray de- 
projection of the cluster temperature and density structure is per- 
formed in section |5] Finally, we incorporate all available SZE 
observations of Abell 2163 in section|6]to derive the best-fit line- 
of-sight peculiar velocity and central Comptonization from the 
SZE spectrum. We list our conclusions in section |7j 

2. Observations 

2.1. APEX-SZ observations at 150 GHz 

Observations of Abell 2163 at 150 GHz were carried out in April 
2007 with the APEX-SZ bolometer camera, with a measured 



beam full-width half-maximum (FWHM) of 58" (determined 
from observations of Mars) and an effective bandwidth of 23 
GHz (Halverson et al. 2008, hereafter H09). At the time of ob- 
servation, 280 of the 330 bolometers were read out, and of these 
about 160 were used for the map making after applying a noise 
cutoff (section[3]l. 

To obtain a uniform coverage over the extended emission of 
the cluster, horizontal raster scans were made with a 15' x 15' 
scan pattern, tracking the central cluster position. The median 
scan speed was 230" s"', and the total integration time was 7.6 
hours. The observing conditions were good to moderate for the 
site, with 150 GHz zenith opacities in the range 0.04-0.09. 

To test a different scan strategy, 0.8 hours of additional obser- 
vations were carried out in August 2007, scanning in circles of 
radius 6 arcminutes. At the beginning of each 2-minute subscan, 
the circle center was placed on the cluster center, and during the 
subscan the source was allowed to drift through the scan pattern. 

2.2. LABOCA observations at 345 GHz 

After producing a high signal-to-noise detection with APEX-SZ, 
five hours of follow-up observations at 345 GHz with LABOCA 
were completed in September 2007. Seven further hours of ob- 
servations at this frequency were carried out in May 2008. 

At the time of observation, LABOCA had 268 optically ac- 
tive bolometer channels, of which around 240 were used for 
mapping of Abell 2163. From observations of Neptune the co- 
added beam size was determined to be 19.5" (FWHM). The ob- 
serving conditions were stable in both observing periods, with 
a typical atmospheric zenith opacity of 0.2-0.3 at the observing 
frequency. Abell 2163 was observed in an elevation range of 40- 
70 degrees. 

For the LABOCA observations, a spiral pattern was used 
with inner and outer radii of 120 and 180 arcseconds, respec- 
tively. Each scan was made up of four such spirals, separated by 
140 arcseconds in azimuth and elevation, to make the coverage 
close to constant over the central part of the cluster 

2.3. X-ray observations 

X-ray emission originating from clusters provides additional in- 
formation on their ICM, and in the interpretation of our results 
we make use of archival XMM observations of Abell 2163. The 
cluster was observed in August 2000 as part of the guaranteed 
time from the first year of operation. Due to the large size of the 
target, as compared to the XMM field of view, a mosaic of five 
pointings was required to fully cover the cluster: one 'on-source' 
pointing^ and four offset pointing^ to probe the outskirts of the 
cluster and the surrounding cosmic background. The nominal 
exposure time of each observation was about 30ks. Since the 
full exposure could not be achieved for the 'on-source' point- 
ing, the observation was completed with a second pointing in 
September 2001, which is included in our analysis. 

3. Data reduction, calibration and mapping 

3.1. General sky-noise considerations 

The SZE signal from Abell 2163 is extended over more than 
20 arcminutes, which is comparable to the field-of-view of 

' Observation Id: 0112230601 

- ObsId: 0112230701, 0112230801, 0112230901, 0112231001 
^ ObsId: 0112231501 
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APEX-SZ (23' across) as well as that of LABOCA (11' across). 
This makes the subtraction of atmospheric signals very difficult; 
subtracting common-mode signals across the array and apply- 
ing polynomial baselines (corresponding to large-scale Fourier 
modes) in this case has the side effect of also removing astro- 
nomical signals on scales comparable to the field-of-view of the 
array. 

The properties of the two instruments used in this anal- 
ysis differ substantially, calling for highly specialized reduc- 
tion schemes in each case. In general, the APEX-SZ observa- 
tions were carried out in poorer observing conditions than the 
LABOCA observations, with relatively high levels of precip- 
itable water vapor. For this reason, the APEX-SZ data suffers 
from excess low-frequency noise correlated on scales smaller 
than the array, requiring high-pass filtering of individual bolome- 
ter time streams to be applied after removing the correlated at- 
mospheric signal. While this step enhances the signal-to-noise 
ratio of the detection, it also removes additional astrophysical 
signal. To account for this, we make use of a point source trans- 
fer function as in H09. It is used as follows: (i) for the isother- 
mal analysis, we model the cluster by convolving the param- 
eterized cluster model with the transfer function and compare 
the result with our reduced map; and (ii) for the non-parametric 
non-isothermal analysis, we deconvolve the reduced map with 
the transfer function. 

In spite of the higher frequency band, there is no excessive 
low-frequency atmospheric noise component in the time streams 
of the LABOCA data. However, because of the lower sensitiv- 
ity of this instrument with respect to the SZE signal, a smaller 
scan pattern was used, compared to the APEX-SZ observations, 
to maximize the signal-to-noise ratio in the limited time avail- 
able. The small scan pattern fundamentally limits the scale on 
which the SZE signal can be recovered. In addition, the details 
of the LABOCA data reduction limits the use of a transfer func- 
tion since the latter cannot be constructed to be linear; in other 
words, the reduction pipeline acts differently on a point source 
and an extended source. To recover as much of the SZE incre- 
ment signal as possible, we use a slightly modified form of the 
iterative map-making algorithm outlined by Enoch et al. (2006). 

For the complete reduction, the bolometer array data analysis 
software BoAQ has been used. 

3.2. Calibration of tiie millimeter data 

Primary flux calibration for APEX-SZ and LABOCA relies on 
observations of planets. The APEX-SZ data is calibrated with 
daily raster scans of Mars on all the 280 elements of the array (cf. 
H09). From these observations the gain, beam shape and position 
on the sky are determined for each bolometer After adjusting 
the bolometer relative calibrations and angular offsets derived 
from the Mars observations, a co-added rms weighted map is 
made using all optically active bolometers. This map is used to 
determine the overall calibration factor, taking into account the 
significant side lobes of the beam. 

The calibration procedure is similar for LABOCA, but here 
the primary calibrator is Neptune, which remains unresolved in 
the 19.5" beam. Because Neptune was not always observable, 
several secondary calibrators were used for LABOCA (Siringo 
et al. 2009). These sources used have been monitored by the 
LABOCA team since early 2007, and their 345 GHz flux densi- 
ties, relative to Neptune, are known to within about 6%. 



The absolute flux of Mars at the APEX-SZ frequency at each 
observing period is determined using a modification of the Rudy 
model (Rudy et al. 1987; Muhleman & Berge 1991), maintained 
by Bryan Butlef] and corrected to be in agreement with recent 
WMAP results as described in detail in H09. The absolute flux 
of Neptune at 345 GHz is derived from cross-calibrations of 
Neptune against Mars (Griffin & Orton 1993). 

To monitor the stability of opacity corrections at different el- 
evations, secondary-calibrator measurements are used. The com- 
putation of opacities relies on skydip measurements combined 
with tau-meter readings (Weiss et al. 2008, Siringo et al. 2009). 

3.3. APEX-SZ data reduction 

Time stream data from APEX-SZ is processed through a data 
reduction pipeline and binned to form maps. This process is de- 
scribed here. The pipeline is slightly different for raster scans 
and circular drift scans, but identical for each scan within these 
two subsets. 

After eliminating detectors with low optical response, corre- 
lated (atmospheric and electronic) noise is removed by subtract- 
ing the median signal from across the good channels of the array 
at each time sample. Individual time streams are then despiked 
by flagging and removing data that deviates from the baseline 
by more than 5cr, which typically corresponds to less than 0.1% 
of the data. Electronic glitches are recognizable because they oc- 
cur with durations shorter than the detector response time. These 
features, seen only in a negligible amount of data, are removed 
as well. 

To baseline the data of a raster scan, the scan is divided 
into subscans of constant elevation, resulting in subscans ex- 
tending over 30' in azimuth. Turnarounds in the scan pattern, 
where the scan speed is significantly lower, are flagged and re- 
moved. A fifth order polynomial is fitted to and subtracted from 
each constant-elevation subscan, effectively high-pass filtering 
the time streams. 

Circular drift scans are baselined by defining subscans con- 
sisting of 2.5 full circles and applying a fifth order polynomial, 
after flagging the first full circle of the scan (which usually con- 
tains intervals of high scan acceleration). This polynomial base- 
line corresponds to a spatial filter similar to that applied to the 
raster scans. Simultaneously, an airmass correction is applied as 
described in H09. 

For each scan, a map with 10" x 10" sized pixels is con- 
structed, weighting the data by the inverse rms of each reduced 
time stream. All maps are then co-added. In parallel with the 
reduction of the data, a beam-shaped source, without noise and 
translated into time stream data, is passed through an identical 
pipeline to obtain the point source transfer function, as described 
by H09. All flags and weights on the data are carried over to the 
artificial data. 

The resulting transfer function is used to deconvolve the high 
significance co-added map of Abell 2163 to the intrinsic reso- 
lution of the instrument. The deconvolution is done iteratively, 
and is similar to the CLEAN algorithm for interferometry data 
(Hogbom 1974, Schwarz 1978) in that the source is modeled as a 
sum of many point sources; here, however, the process is carried 
out in map space rather than Fourier space. The raw (pipeline 
filtered) map, Mq, is first convolved with the beam to reduce the 
sample variance on small spatial scales. Each pixel in the map 
is then divided by the local rms to create a signal-to-noise map, 
A^o, from which the brightest pixel is selected. The correspond- 



'* http://www.apex-telescope.org/bolometer/laboca/boa/ 



' http://www.aoc.nrao.edu/ bbutler/work/mars/model 
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ing point source flux corresponding to this peak pixel value is 
computed, taking into account both the beam smoothing of Mq 
and the complete description of the inherent resolution and the 
effects of the reduction described by the point source transfer 
function. The resulting point source is re-convolved with the 
transfer function to represent how the point source is reflected 
in the raw map. This component is subtracted from Mq to yield 
Ml, from which a new smoothed signal-to-noise map A^i is con- 
structed. The point source is also convolved with the "clean" 
beam and added to a map C of "clean component". Next, the 
brightest pixel in A^i is located, and the process is repeated un- 
til Mi is consistent with noise. The criterion for the latter is that 
as many negative as positive pixels are selected in any consec- 
utive 20 iterations. At this point, M, is added to the cumulative 
map C to account for emission which has not been successfully 
removed as point sources. Details on the limitations of this de- 
convolution technique will be given in a future publication (Nord 
et al., in preparation). The de-convolved APEX-SZ map of Abell 
2163 is shown in Figure [1] 

To compute the noise on the beam scale, we smooth the fi- 
nal map with the APEX-SZ beam (FWHM 58" and computing 
the rms in square regions with side five times the beam. In this 
way, we find a noise rms value of 32 fiKcuB on the beam scale, 
corresponding to 0.013 MJy sr in the central region of the de- 
convolved map. 



from this difl^erence, tiie time stream reduction sequence used in 
each iteration is quite similar to that used for the APEX-SZ data, 
and basically follows the method outlined by Weiss et al. (2008) 
and detailed by Siringo et al. (2009). Maps are constructed with 
6" X 6" pixels, and co-added in the same way as for the APEX- 
SZ data. 



To account for sky signal attenuation in the time stream re- 
duction, we follow the iterative approach of Enoch et al. (2006), 
with one important modification. Using a suitably chosen time 
stream reduction algorithm, a map is produced. Each pixel with 
a significance of less than 3 sigma is set to zero, and the rest of 
the map is chosen as the template to be subtracted directly from 
the time streams before running the same time stream reduction 
yet again, and so on. While Enoch et al. take a conservative ap- 
proach and derive the next template from the residual signal, we 
add the template back into the data prior to mapping and derive 
this template from the total map. This ensures that any signifi- 
cant feature in the current best guess of the source flux distribu- 
tion is carried over to the next iteration step. 



SZE decrement (nK CMS) 
-0.6 -0.4 -0.2 0.0 
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Fig.l. Map of Abell 2163 at 150 GHz, overlaid with 
XMM-Newton X-ray contours (see Fig. |3]l in units of 
10^'^erg s^'cm^^arcmin"^. Because the correlated-noise re- 
moval has attenuated the source signal, the map has been de- 
convolved using the point source transfer function (see text). 



3.4. LABOCA data reduction 

Because the LABOCA time streams are more stable at low fre- 
quencies than those from APEX-SZ, the data do not have to 
be further filtered after standard correlated-noise removal. Aside 



The iterative mapping algorithm requires 8 iterations before 
convergence is reached. In each iteration, an identical time series 
reduction is performed. 



Because the iterative mapping technique is non-linear, it is 
not straightforward to characterize it by a filter function. In place 
of a transfer function, the instrument beam is used when fitting a 
model to the 345 GHz data. Significant loss of signal is expected 
using the iterative approach. To quantify the level of bias, we use 
the best-fit fi model from the APEX-SZ measurement (section 
and pass is through the complete LABOCA data reduction 
pipeline. Although more than 30% of the signal is lost beyond 
'"500 (the radius at which the enclosed average matter overden- 
sity is 500 times the mean cosmic density), we find that the sig- 
nal loss within r25oo (defined analogously, approximately 3') is 
negligible compared to statistical errors. The central rms in the 
LABOCA map, computed on a scale of five times the LABOCA 
FWHM (19.5"), is 1.9 mJy/beam, corresponding to 0.24 MJy 
sr-i. 



We find one significant point source in the final map (Figure 
|2|l, likely thermal emission from a cluster galaxy or a background 
source lensed by the cluster. After highpass-filtering the map 
to remove structures larger than two beams (such as the cluster 
emission), we find a peak flux density of 11.8 ± 1.9 mJy/beam 
for this source. 
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Fig. 2. Top: Final 345 GHz LABOCA map of Abell 2163, 
smoothed with the 19.5" beam. The cross marks the position 
of a bright flat-spectrum radio source (Cooray et al. 1998). The 
diamonds mark the positions of two BCGs (Maurogordato et al. 
2008). The circle marks the position of the bright point source 
found in the present data. Bottom: LABOCA map smoothed to 
the APEX-SZ resolution of 1 arcminute. The bright point source 
described in the text has been removed. The APEX-SZ 150 GHz 
map is shown as contours. 



3.5. XMM-Newton data reduction 

Here we describe flie processing steps applied to the X-ray data. 

We generate calibrated event-lists from the raw data of the 
6 archival pointings using the standard procedures of the XMM 
Science Analysis Systen^ (SAS) and screen these lists for pe- 
riods of high radiation due to soft solar proton flares. To this 
end, we apply the method of Pratt & Arnaud (2002), which con- 
sists of fitting histograms of the high energjQ light curves with 
a Poisson law and rejecting time intervals exceeding the mean 
radiation level by more than 3cr. On average, this shortens the 
available exposure time by 15%. Although the high energy sig- 
nal is very sensitive to particle flares, one needs to adopt a coarse 
time sampling (104s here) to ensure that the average number of 
counts per bin is significant. This implies that the wings of light 



curve jumps in increased radiation periods cannot easily be de- 
tected. To improve on our first filtering, we thus repeat the same 
analysis on a broad low energy band, namely [0.3-10] keV with 
narrow time bins of 13s, resulting in the exclusion of an extra 3% 
of the data. Out-of-time event lists are also generated using the 
simulation mode of the SAS tasks emchain/epchain and filtered 
using the same time intervals. 

From the filtered event lists, images and exposure maps are 
generated in the [0.5-2] keV band and gathered into two large 
mosaics. The signal to noise of the cluster in this band is about 
80% of the best achievable value using the XMM bands. It 
is however free from the strongly varying instrumental lines 
around 6-7 keV and calibration uncertainties at lower energies. 
Combined with the weak dependence of the cooling function on 
temperature at low energies (soft X-rays), this makes it a very 
suitable band for our purpose. Since the XMM response is al- 
most flat in this narrow band, the exposure maps are evaluated at 
a single average energy of 1 .25 keV. 

The main issue in the X-ray analysis of extended source is 
an accurate background modeling. To this end, we make use of 
the Filter Wheel Closed (FWC) dataset provided by the XMM 
background working groupQ. Our model consists of a sum of a 
re-scaled FWC image (which represents the instrumental noise), 
an homogeneous cosmic background of constant surface bright- 
ness (since the area of the mosaic remains small) and a properly 
scaled image of the out-of-field events. The scaling of the instru- 
mental background is allowed to vary with each pointing and 
instrument while the cosmic background level is just a function 
of the detector (to account for the different responses). Using the 
Cash statistics (Cash 1979), the model parameters are fitted si- 
multaneously to the 6 observations including the data from the 
out of field corners and excluding the positions closer than 17' 
from the cluster center. 

In order to detect and mask the surrounding AGNs, we adap- 
tively smooth the background subtracted raw mosaic, using the 
criterion of having more than 10 counts per cell. The resulting 
map is then exposure corrected and smoothed again by a cr = 2" 
GaussiarQ A source catalog is extracted from this image using 
SExtractor (Bertin & Arnouts 1996) and the associated source- 
segmentation map is combined with an exposure time threshold 
in order to derive a global mask. 

The tasks of background estimation and source detection are 
somewhat intricate, especially in the case of mosaics where steep 
instrumental and particle background variations can occur in be- 
tween pointings. To tackle this issue, we iteratively re-perform 
the background modeling and AGN extraction using the updated 
AGN mask for the background estimate. Both products stabilize 
after three iterations. 

The final background subtracted, adaptively smoothed image 
of the cluster is shown in Fig. [3] 



4. Isothermal modeling of the intra-cluster gas 

In order to assess our SZE measurements quantitatively and 
compare them with the X-ray data, we fit the SZE measure- 
ments using the well-known isothermal /3 model (Cavaliere & 
Fusco-Femiano 1978). This enables us to compare the X-ray 



* http://xmm.esac.esa.int/sas/, all the data processing for this paper 
relies on v7.1.2 

^ more precisely [10-12] keV for the MOS detectors and [12-14] keV 
for the PN 



http://xmm2.esac.esa.int/external/xmm_sw_cal/background/ 
' This hybrid smoothing allows to wash out the Poisson fluctuations 
in regions devoid of source signal without degrading the resolution over 
moderately bright sources (which are both the advantage and drawback 
of adaptive filtering alone) 
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Fig, 3. Smoothed, background subtracted. X-ray map of 
Abell 2163 in the [0.5-2] keV band (see text for details). 
Logarithmically spaced contours highlight the broad dynami- 
cal range of the cluster emission. The unit of the color scale is 
erg s^'cm^^arcmin"^. 



surface brightness profiles with the SZE temperature decre- 
ment/increment maps and check the agreement of our SZE mea- 
surements with the prediction from the X-ray mean values. 

4.1. Isothermal gas in hydrostatic equilibrium 

The isothermal /3 model is widely used to model the gas profile 
in clusters under the assumption of hydrostatic equilibrium. The 
radial distribution of gas density is given by 



(1) 



n,(r) = n,(0) 1 + - 



where is the density of electrons, r is the radius from the cen- 
ter, is the core radius of the gas, and /? is the power law index. 

With this model, the temperature decrement (or increment) 
due to the SZE takes the form 



AT(e) = fix, T,) TcMB y = AT(0 = 0) 1 -H 



,2\(1-W2 



0? 



(2) 



where 6 is the angular radius on the sky, AT" is the thermody- 
namic SZE temperature decrement (or increment) and y is the 
Comptonization parameter, given by the line-of-sight integral 



y 



J [meC^j 



CTj dl. 



(3) 



where crj is the cross section for Thomson scattering. The fre- 
quency dependence f(x, Te) of the SZE signal is given by 



+ 1 



e" - 1 



(1 +<5.,r,), 



with relativistic corrections contained in 6xj, as described by 
Sazonov & Sunyaev (1998). x is the dimensionless frequency 



The analogous expression for the X-ray surface brightness 
reads 

A5x(0) = 5„ 1 + -^ , (4) 



where 5o is the central surface brightness and the different ex- 
ponent comes from the dependence of X-ray emissivity. 

4.2. Elliptical p model fit to the SZE data 

Because of the significant ellipticity in the X-ray surface bright- 
ness profile, we generalize equation|2]to an elliptical form given, 
e.g., in H09 as 



ArszE = Aro(l+A + B) 



(l-3/?)/2 



(5) 



where 



A = 



(cosmiX - Xq) + sinmiY - Yo)f 



B 



(- sin((I>)(X - Xq) + cos((D)(7 - Ya)Y 



Here {X - Xq) and (Y - Yq) are angular offsets on the sky in the 
right ascension (RA) and declination (DEC) directions, respec- 
tively, with respect to centroid position (Xo, Yo). The axial ratio, 
T], is the ratio between the minor and major axis core radii and <1> 
is the angle between the major axis and the RA (X) direction. 

An elliptical /3 profile is fitted to the 150 GHz data by con- 
volving each model function with the point source transfer func- 
tion and minimizing ax^ statistic weighted by the inverse square 
of the local map rms. Note that the raw, not the de-convolved, 
map is used for this fit in conjunction with the transfer function. 
All map pixels within a 10' radius of the X-ray center (section 
I4.3l l are considered for the fit. The results of the fit are given in 
Table □ 

For comparison, we also fit a spherical /5 profile and find 
roughly consistent results. At 345 GHz, the detection of Abell 
2163 is less significant, we thus fix the centroid and the value of 
/? to those fitted from the 150 GHz measurement and fit for the 
remaining parameters. 

Errors in the fitted profiles are estimated using a Monte Carlo 
approach. The map data is jack-knifed by inverting the sign of 
half the individual scan maps (randomly chosen), and the maps 
are co-added to generate a pure noise map. After convolving the 
best-fit model with the transfer function and adding the result 
to the noise map, the statistic is again minimized to find the 
best set of parameters for the new map. The entire procedure 
is repeated 1000 times, upon which the errors in the respective 
parameters are taken as the scatter in the distributions of their 
fitted values. We verify that the mean of the distributions of fitted 
parameters are consistent with the best-fit model. 

4.3. Elliptical p-model fit to the X-ray data 

To allow for a basic comparison of the SZE and X-ray data in the 
context of the isothermal /3 model, we perform an isothermal fit 
to the X-ray data, similar to that used for the SZE data. As for the 
SZE analysis, we generalize the spherical /3 model of equation|4] 
to the elliptical case, as: 



ASx = A5o(l +A + B) 



1/2-3/; 



(6) 
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This model is fitted to the raw X-ray images using the Cash 
statistic (Cash 1979). The Poisson nature of the data does not al- 
low here for a direct background subtraction; instead we add the 
simple two-component background described in section 13.51 to 
the fitted components. The results of the fit are presented in Table 

4.4. Systematic uncertainties 
4.4.1. Galactic dust 

The effect of galactic dust emission on the SZE measurement 
of Abell 2163 has been estimated by Lamarre et al. (1998) for 
the PRONAOS balloon experiment and by LaRoque et al. (2002) 
for the SuZIE experiment (Holzapfel et al. (1997a,b)). However, 
because these two estimates apply to scan strategies and spatial 
filtering functions very different from ours, it is not appropriate 
to extrapolate from either of them. 

To estimate the dust correction at 345 and 150 GHz, we use 
the IRAS far-infrared 60 //m and 100 jim dust maps (Schlegel et 
al. 1998) on the Abell 2163 field . Smoothing the 60fim map to 
the 4.3' resolution of the 100 /^m map and re-sampling the maps 
so that each pixel has a side of one arcminute, we extrapolate 
to the frequencies of the SZE measurements using a grey-body 
spectrum of dust emissivity 



exp(hv/kBTd) - 1 ' 

where a is the dust spectral index. Using a - 2.0 (Finkbeiner 
et al. 1999) we fit to for every 1' x 1' pixel pair in the two 
dust maps. After smoothing the resulting temperature map with 
a 4.3' Gaussian kernel to remove noise artifacts, we use equation 
(|7|l to extrapolate the intensity at 345 and 150 GHz, normalizing 
at lOOyum. 

The simplest way of removing the dust component from the 
data is to subtract it before any filtering is applied by the re- 
duction. We thus create time streams corresponding to the ex- 
trapolated dust maps on the coordinates covered by APEX-SZ 
and LABOCA and subtract these from the raw time streams be- 
fore proceeding with reductions identical to those described in 
section 13 In this way, we find corrections of -0.6% to the SZE 
decrement at 150 GHz and -1-1.8% to the SZE increment at 345 
GHz. Effects on the other fitted parameters are negligible. We 
assume that contamination on angular scales smaller than four 
arcminutes (i.e. not resolved by IRAS) are also negligible. 

4.4.2. Point source contamination 

After subtracting a preliminary (possibly point source contami- 
nated) /3 model from each map, we smooth the maps to the re- 
spective beam scales and look for point source detections above 
3cr. One source (cf. Figure |2]i is found in the LABOCA map at 
345 GHz. No sources are found in the APEX-SZ 150 GHz map. 
Point sources which are detected with sufficient S/N in the maps 
could be removed by subtracting fitted Gaussians. Here we take 
a somewhat more conservative approach and simply disregard a 
region around the 345 GHz point source in a subsequent fit with 
a new /3 model. 

Extragalactic infrared point sources are not expected to cor- 
relate spatially with clusters of galaxies, and a distribution of 
sources can be well approximated by excess noise in any dif- 
ferential measurement (White & Majumdar 2004), thus having 
little systematic effect on bolometric SZE measurements. 



Although tensing by a massive cluster such as Abell 2163 
can significantly increase the number of detected sources in the 
cluster central region, this effect conserves the total intensity of 
the background and only becomes important when faint sources 
are raised above the detection limit and removed, and the un- 
resolved background intensity is thereby systematically lowered 
(e.g. Loeb & Refregier 1997). It is not clear whether the bright 
point source found in the 345 GHz map is a tensed background 
source; however, we have verified that the effect in fitted pa- 
rameters of excluding the region around this source for the fit 
is merely 4%, i.e. smaller than both the derived statistical errors 
and the absolute calibration uncertainty. 

Radio sources are expected to correlate strongly with clus- 
ters of galaxies (e.g. Reddy & Yun 2004). We have searched the 
NVSS catalog (Condon et al. 1998) for sources brighter than 
5 mJy at 1.4 GHz. The flux of each source is extrapolated to 
150 GHz using the spectral index maps of Feretti et al. (2004). 
With this simple approach, even the brightest source will have a 
peak flux considerably lower than the rms in the APEX-SZ map. 
Cooray et al. (1998) report an inverted-spectrum non-thermal 
source of 3 mJy at 30 GHz, indicated in Figure|2] No counterpart 
can be seen in either SZE map. High-pass filtering the APEX-SZ 
map, after subtracting the best-fit model of the cluster, to remove 
all structure larger than 3' results in an effective rms of 30 yuKcMB 
at the position of this source. The fact that the source is not seen 
in this filtered map allows us to put a lower limit on the spectral 
index as qt^q" > 0.5 (5y oc v^"). Subtracting the corresponding 
signal from the APEX-SZ map before the /? model fit results in 
a systematic shift in the central SZE decrement of less than 1%. 
Based on these considerations, we consider systematic effects 
from radio sources negligible in comparison to primary CMB 
anisotropics and absolute calibration. 

4.4.3. Primary CIVIB contamination 

For a very extended cluster Uke Abell 2163, with an estimated 
rsoo greater than 7' (Radovich et al. 2008), the temperature 
anisotropics in the CMB are a potential source of confusion in 
the SZE measurements. This uncertainty particularly affects the 
kinematic SZE signal because the latter has a frequency depen- 
dence identical to that of the CMB thermal spectrum. Note that 
this is not an issue for CMB power on scales > 15' due to the in- 
direct high-pass filtering of our maps (through polynomial base- 
lines in the bolometer time series). 

To quantify the level of contamination on scales compa- 
rable to the SZE emission, the HEALpix software (Gorski et 
al. 2005) is used to generate 100 realizations of the CMB sky 
in a one square degree field with a resolution of 1.7'. The in- 
put power spectrum for this generation is computed using the 
CMBfast code (Zaldarriaga & Seljak 2000) with the WMAP 5- 
year cosmology (Komatsu et al. 2009). Relying on the results 
of LaRoque et al. (2002), diffuse secondary CMB anisotropics 
other than the SZE are considered negligible. 

Each primary-CMB map is re-sampled as time streams cor- 
responding to a typical scan of the SZE observations at the two 
frequencies and run through the reduction pipelines to account 
for filtering. The /3 model fit is repeated after subtracting each of 
the thus filtered primary-CMB realizations. The scatter in the re- 
sults of each parameter is taken as the 68% systematic marginal 
error from primary CMB contamination. We find a fractional er- 
ror of 3.4% at 150 GHz and 2.9% at 345 GHz. 

The primary-CMB component has been simulated on the 
relevant scales compared to the cluster emission, with scales 
smaller than 1 .7' not taken into account. This is justified because 
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the power of the primary anisotropies drops rapidly on scales 
smaller than a few arcminutes. 



4.4.4. Pipeline filter function 

To test how well the point source transfer function represents the 
filtering of the SZE signal in map space at 150 GHz, we carry out 
a series of reductions of known signals. To this end, we remove 
astrophysical signals from the data by dividing the scans into two 
minute segments and randomly inverting half of these. Before 
passing the noise data through the pipeline, a beam-convolved 
y6 model is added to the time streams. The point source transfer 
function is used to reconstruct the y6 model after reduction, and 
the procedure is repeated for a large number of /? models in the 
range 50" < 9^ < 150", 0.5 < yS < 1.4. We find no evidence 
that the level of noise affects the result other than to increase the 
scatter in fitted parameters; we thus use bright model sources 
(~ 10 mKcMB at 150 GHz) to estimate the effect of the pipeline. 

Although the fitted central amplitude can be systematically 
reduced by as much as 12% for combinations of small /? and 
large Oc (characterizing very broad profiles, where some modes 
in the map are irrecoverable by the transfer function), we find the 
systematic shift to be no more than 2% for either of the models 
fitted to the actual 150 GHz data. Due to the degeneracy between 
y6 and 9c, there is considerable scatter in these parameters when 
fitted to the simulated data. Keeping one parameter fixed and fit- 
ting the other, we find that the systematic errors in both param- 
eters due to inaccuracies in modeling the pipeline filter function 
are negligible. 

4.5. Results of the isothermal modeling 

Figure |4] shows the radial profile of the best-fit spherical /3 pro- 
file to the 150 GHz data and the same profile convolved with 
the transfer function, as well as the profiles of the raw and 
de-convolved maps. The best-fit /? model parameters, corrected 
for systematic effects as described below, are given in Table [1] 
Although derived from a completely independent method, the ra- 
dial profile of the de-convolved map is fully consistent with the 
beam-smoothed best-fit model within the 10' truncation radius 
used for the fit. 

Comparing the results of the X-ray and SZE fits, a non- 
negligible offset (21 ± 8"), possibly caused by an asymmetric 
temperature distribution near the center, appears between the 
two centroid positions. The origin of this offset will be discussed 
as part of further analysis in a future publication (Kneissl et al., 
in preparation). The other fitted parameters show an overall con- 
sistency in the global shape of the emission. We note that in spite 
of the significant ellipticity in the X-ray surface brightness, the 
elliptical model of the SZE signal is only barely inconsistent 
with the spherical one (the axial ratio deviates from 1 by only 
1.7cr). The most likely explanations for this are the poorer reso- 
lution of the SZE data which results in a greater uncertainty on 
the ellipticity, and the inherent fact that the SZE signal, due to 
the density weighting, is naturally smoother and more diffuse. 



5. Non-isothermal modeling of the intra-cluster gas 

Several authors have discussed the joint modeling of SZE and X- 
ray data for de-projecting ICM parameters (see, e.g., Yoshikawa 
& Suto 1999, Zaroubi et al. 2001, Puchwein & Bartelmann 2007, 
Ameglio et al. 2007). However, most of these analyses have 
been limited to analytic cluster models or numerical simulations. 




4 6 8 10 

Angjlar distance (orcminutes) 

Fig. 4. Radial SZE profile of Abell 2163 at 150 GHz. EiTor bai's 
indicate the profile computed from the reduction-attenuated raw 
map (before deconvolution; this is the map used for the para- 
metric fit) while the shaded region represents the profile from the 
map de-convolved to the beam resolution (note that this map was 
not used to derive the best fit /3 model. The profile of the best-fit 
isothermal /3 model is indicated by solid lines; the thick line rep- 
resents that model convolved with the beam, while the thin line 
represents the same model convolved with the transfer function. 
The dashed line indicates the profile of the de-convolved map, 
directly re-convolved with the transfer function. The vertical ar- 
rows indicates the radial cut of 10' used for the model fit, as well 
as r2oo and rsoo. 



Our work here represents the first attempt at de-projecting clus- 
ter density and temperature profiles using actual SZE imaging 
data, in combination with the X-ray data, without resorting to 
any parametric models for the gas or dark matter distribution. 

Given the large angular size of Abell 2163, the arcminute 
resolution of APEX-SZ is sufficient to carry out a joint radial 
density and temperature modeling of the ICM in combination 
with publicly available X-ray data. Although the LABOCA map 
has more than a factor of two better resolution, due to its high 
noise level and small field of view we restrict the analysis in this 
section to the APEX-SZ data. 



5.1. Method 

5.1.1. Deprojection method 

For the deprojection analysis we use a direct inversion technique 
based on the Abel integral. The technique was first proposed by 
Silk & White (1978; see also Yoshikawa & Suto 1999). Apart 
from the assumption of spherical symmetry, there are no ad- 
ditional theoretical constraints (e.g. hydrostatic equilibrium or 
polytropic gas index), making this is a natural method for ex- 
tracting ICM parameters, assuming that spherical symmetry in- 
deed provides a good estimate of the structure. The details of this 
method and its limitations will be discussed in a future publica- 
tion (Basu et al., in preparation); here we present the outline of 
the method and the results for Abell 2163. X-ray spectroscopic 
measurements have shown that Abell 2163 has a complicated 
temperature structure at the center and is most likely a merg- 
ing system (Elbaz et al. 1995, Govoni et al. 2004). Nevertheless, 
the de-projection analysis with spherical symmetry implemented 
here shows how resolved SZE images of clusters can immedi- 
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Table model fit results at 150 and 345 GHz. 



Parameter 


ljU Ljrlz (^eilipticaij 


ijU Oriz ^spnencaij 


j4j Oriz (^spnericaij 


X-ray (elliptical) 


Xo (Central R.A. [J2000]) 


16''15"'45.P±9" 


16M5"45.6»±8" 


(16M5™45.6»)* 


16''15"' 46.6^+1.5" 


Yo (Central DEC. [J2000]) 


-06°08'31"±8" 


-06°08'28"±8" 


(-06°08'28")'' 


-06°08'44.1"+1.4" 


Peak signal" 


-0.301+0.033 


-0.302+0.033 


0.613+0.089 


1.075 X 10-'2+7 X 10^'' 


9c (core radius) 


123±19" 


129+21" 


89+29" 


90.7+2.8" 


P (power-law index) 


0.85+0.14 


0.85+0.14 


(0.85)* 


0.639+0.03 


<t> (inclination angle) 


128±24° 


(0.0)* 


(0.0)* 


1.2+6° 


ri (axial ratio) 


0.77+0.11 


(1.0)* 


(1.0)* 


0.850+0.05 



" Brightness in units of MJy sr ' (SZE), and in units of erg s ' cm - arcmin - (X-ray) 
* Fixed parameters indicated by parentheses 



ately be used in combination with X-ray data to gain a better 
understanding of the gas and mass distribution. 

The deprojection analysis is based entirely on the radial pro- 
files of X-ray surface brightness and SZE temperature decrement 
(or increment). From equation ^ the SZE temperature decre- 
ment can be written as the integral over the line of sight; 



Jr 



AT(R) = 2Asz fix, Te) iieir) Te(r) 



rdr 



Vr2 - 



(8) 



where the cluster is assumed to be spherically symmetric, Asz - 
(Tj (kTcMB hn-eC^), r is the physical radius from the cluster center, 
R - Da9 where 9 is the projected angular distance on the sky 
and Da is the angular diameter distance (to calculate Da, we 
assume the cosmology of Komatsu et al. 2009 for the remainder 
of the paper). Te{r) and ne{r) are the electron gas temperature 
and density radial profiles. We neglect the small dependence 
in f(x, Te) for this analysis and incorporate f{x, Te-\0 keV) into 
the Asz factor 

For the de-projection analysis, the X-ray surface brightness 
profile can be written as (e.g. Yoshikawa & Suto 1999): 



SxiR) 



(r) KniTAr)) 



rdr 



Vr2^ 



R-^ 



(9) 



We use the APEC code (Smith et al. 2001) to compute the 
cooling function, MJeif)) for Abell 2163, assuming an abun- 
dance value of 0.4Zq, where Zq is the solar metallicity. This also 
includes the absorption by the neutral hydrogen of the Galaxy 
using a column density of 1.1 x 10^' cm"^, as measured by 
the L.A.B. survey (Kalberla et al. 2005) at the cluster position. 
While the temperature dependence of the cooling function is 
thus included in the analysis, it is interesting to note that the 
results changes very little when assuming A is independent of 
Tg. As we shall see in section |5^ however, systematic effects 
in the gas density from uncertainties in the metal abundance are 
non-negligible. 

Using Abel's integral equation, equations ([8]) and (|9]l can be 
inverted to obtain joint radial density and temperature profiles: 

1 f dATjR) AR 
Te{r)ne{r)^— — — — p===; (10) 

%J oo 



;rA 



sz 



AR 



n^(r)A(r,(r)) = 4(l+zy 



\J CO 



ASxiR) dR 



dR v^2~r72 



(11) 



Following Yoshikawa & Suto (1999), we integrate equations 
( [Tol l and (fTTT i numerically by summing in radial bins from /min to 
'max, where /^ax is the index for the outermost bin, and corre- 
sponds to r/DA- We propagate the errors on the density and tem- 
perature profiles by a Monte-Carlo method, the details of which 
will be discussed in a future publication (Basu et al. 2009, in 
preparation). 



5.1.2. Profile extraction 

The SZE profile is derived as a radial average of the data and 
errors are estimated using the same method as in H09. In order 
to combine the two data sets, the X-ray raw mosaic has been 
smoothed to the APEX-SZ resolution of 1'. In doing so, we have 
neglected the effects of the XMM point spread function (PSF) 
and its variation over the field since they are negligible com- 
pared to the APEX-SZ beam. The X-ray source mask is used to 
remove the sub-cluster 8 arcminutes north of Abell 2163, promi- 
nent in Figure|3] and the numerous AGN. A count profile is then 
extracted from the background subtracted mosaic and corrected 
for the local exposure time and unmasked area, thus converting 
to surface brightness. The Poisson error bars on the counts pro- 
files are re-scaled accordingly. Since the background and source 
counts in each annulus are large, the noise can be considered 
as Gaussian which justifies the direct background subtraction. 
Because of the slight offset between the peaks of the X-ray sur- 
face brightness and the SZE temperature decrement, we have 
taken the centroid position of the X-ray map to be the common 
center in the X-ray/SZE joint analysis because of the much better 
signal-to-noise ratio of the X-ray map. 



5.1.3. Mass estimation metliod 

Using the SZE and X-ray measurements, we can determine the 
gas mass and total mass enclosed within a certain radius assum- 
ing hydrostatic equilibrium. This can be used to compute the 
gas-to-mass ratio as a function of radial distance, or can be com- 
bined with the weak-lensing data for a more direct determination 
of the gas mass fraction. We show the comparison between the 
gas mass fraction values obtained using our non-parametric de- 
projection method and the standard isothermal modeling inside 
the rsoo of the cluster. 

Computing the gas mass and total mass profiles from isother- 
mal /3 models is straightforward (e.g. LaRoque et al. 2006). The 
gas mass is obtained directly from the electron density profile as 



Mgas(< r) = Att iJ,ne(0)m^ 



Jo 



el 



eMe, (12) 



where jj^ is the mean molecular weight per electron, which we 
assume to be equal to 1.17 for cosmic abundance of H and 
He. For the non-isothermal analysis, we use the de-projected 
electron density profile to compute the gas mass directly as 

The total mass, Mtotah is obtained by solving the hydrostatic 
equilibrium (HSE) equation, assuming spherical symmetry, as 
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follows: 

ksTeir) r \ dinner) d\nT,(r)] 

Mtotai(< r) = — — + — . (13) 

Ofirrip [ amr amr 

For isothermal modeling, we use the simple analytic form ob- 
tained from the above equation (e.g. Grego et al. 2001) 

M.U<.)=^V^, (14) 

whereas for non-isothermal modeling we solve equation [13] di- 
rectly. The gas mass fraction is in both cases computed as 

/gas(< r) = Mgas(< r)/M,otal(< r). 

5.2. Systematics 

The dominating source of error in the non-parametric de- 
projection analysis comes from uncertainties in the SZE map 
at 150 GHz. Intrinsic noise properties in the raw map cause the 
uncertainties in the radial bins to be strongly correlated due to 
noise structures much more extended than the bin width. In addi- 
tion, the deconvolution process discussed in section [33] not only 
recovers source signal, but can also amplify noise structures in 
the raw map, thus increasing the uncertainties in the radial pro- 
file (see Fig. nil. This amplification of noise eflFectively limits the 
range of radii for which we can make meaningful estimates of 
the ICM temperature and density. 

As already discussed in section 14.4.41 some modes corre- 
sponding to large angular scales are irrecoverable by the transfer 
function - in particular, no structures more extended than the 
path traced by a single bolometer during a scan can be recov- 
ered. Although the SZE profile is well constrained within rsoo, 
this inherent filtering can lead to systematic shifts in temperature 
and density at larger radii, where the SZE signal is potentially 
underestimated. 

Cavaliere, Lapi and Rephaeli (2005) showed that the Abel 
deprojection method in itself can lead to biases in the estimated 
temperature and density profiles. However, they also found that 
this bias is significantly smaller than the statistical uncertainty, 
which they optimistically assumed to be 1 % on the central value. 

Other systematic uncertainties come from primary CMB sig- 
nals, modeling of the X-ray background, and the metal abun- 
dance used in computing the X-ray cooling function. We con- 
sider each source of systematic uncertainty in what follows. The 
results are summarized in table |2] Representative radial scales 
are taken from the best-fit NEW model from weak-lensing anal- 
ysis of this cluster (Radovich et al. 2008): r25oo - 3.1' and 
rsoo = 7.6' . 

5.2.1. The filter function 

To address the issue of how irrecoverable signals on large angu- 
lar scales in the SZE map affect the temperature and density pro- 
files, we once again carry out a series of simulations (cf. section 
14.4.4b . passing artificial models through the reduction pipeline. 
While in section l4~4l we considered only how the signal was af- 
fected by the pipeline in terms of fitting a parametric function 
out to a certain radius using the transfer function, we now also 
have to consider systematic effects from using the transfer func- 
tion to deconvolve the map rather than to convolve a model (the 
latter is much more straightforward). Lacking a more complete 
model of the cluster emission, we again use spherical j6 mod- 
els, and vary the best-fit parameters from section |4] within their 



Icr errors to test for stability. After passing the models through 
the pipeline, we deconvolve them using the transfer function and 
compare the resulting profiles with the input models. While the 
systematic errors estimated from this method are comparable to 
the random errors on the SZE profile inside rsoo, the signal is 
systematically lowered by as much as 40% of the input signal at 
'"200- Note, however, that at any radius considered in our analysis, 
this systematic effect is small compared to the intrinsic statistical 
error originating in the raw map. 

It should be stressed that these error estimates are only mean- 
ingful in the context of the isothermal model. They are justi- 
fied, however, by the results of section|4] which indicate that this 
model indeed provides a very reasonable fit to the data out to 
rgoo. In contrast to what was done in the isothermal analysis, we 
do not attempt to correct for this effect since we have no para- 
metric model. Instead, we treat the systematic as an added uncer- 
tainty on the SZE profile, and propagate this uncertainty through 
the de-projection analysis to estimate the resulting uncertainties 
in temperature and density. 



5.2.2. Primary CMB 

The systematic uncertainty from the primary CMB signal is esti- 
mated using a method analogous to that described in section |4!4l 
After a set of simulated primary-CMB maps are subtracted from 
time stream data and passed through the pipeline reduction, de- 
convolution is performed on each map in the set. The different 
maps are then passed through the de-projection analysis to es- 
timate the impact on the derived quantities. It is found that the 
uncertainties are comparable with those of the filter function. 



5.2.3. X-ray background and metal abundance 

Although the statistical errors on the recovered temperature and 
density profiles are clearly dominated by the uncertainty in the 
SZE signal, several X-ray systematics have to be considered. 

As mentioned in section 13.51 the background model can 
have a significant impact on the X-ray surface brightness at 
large radii. To quantify this further, we fix all the parameters 
of the background model to their best fit value, and allow for an 
overall change in normalization in the instrumental and cosmic 
background. Identifying the Cash statistic with a distribution 
(valid in the large number limit), we derive a Icr uncertainty for 
this normalization factor of -1-1.9%/- 1.2%. At rgoo, tfiis translates 
into a 3% uncertainty in the temperature profile and 1% for the 
electron density, much below the statistical errors. 

Our assumption of a fixed metal abundance is another poten- 
tial source of systematics. Line emission is indeed preeminent in 
the soft energy band and this could result in a significantly differ- 
ent value of the cooling function required for our de-projection 
analysis. We have thus investigated the impact of a global change 
in metallicity in the range 0.2 to O.SZ©. These values correspond 
to the most extreme 1 cr bounds measured anywhere in the cluster 
by Snowden et al. (2008) from X-ray spectroscopy, yet they only 
change our results by up to few per cent, again smaller than the 
statistical uncertainties. The small systematic uncertainties are 
a direct consequence of the high temperatures measured for this 
cluster; line emission at soft energies is indeed much reduced for 
plasma hotter than ~2 keV in which more energetic transitions 
are favored. 
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Table 2. Systematic effects in the de-projected temperature and Angular distance (arcminuies) 

density profiles. Statistical errors are indicated for comparison. o 01 
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Source of 


Effect 


on Tf 


Effect 


on He 


uncertainty 


''2500 


^■500 


''2500 


^500 


Primary CMB 


±3.7% 


±8.7% 


±1.1% 


±3.3% 


Filter function (SZE map) 


+2.1% 


+8.9% 


+0.3% 


+ 1.4% 


X-ray background 


+0.5% 
-0.3% 


+24% 
-3.1% 


+0.8% 
-0.2% 


+1.1% 
-0.3% 


Metal abundance 


+2.2% 
-1.3% 


+ 1.5% 
-0.2% 


+1.9% 
-2.8% 


+2.4% 
-3.5% 


Total systematic_^ 


+4.8% 
-3.9% 


+12.8% 
-9.2% 


+2.4% 
-3.0% 


+4.5% 
-4.8% 


Statistical 


±19.3% 


±43.5% 


±1.1% 


±6.0% 



" Added in quadrature 

5.3. Results of the non-isothermal modeling 
5.3.1. Density and temperature profiles 

The radial profiles of the ICM density and temperature, obtained 
from equations [TO] and (TT] are shown in Figure |5] The density 
profile shows very little deviation from an isothermal /3 model in- 
side rsoo (roughly 1500 kpc), which is expected since the density 
values are mostly constrained by the X-ray surface brightness 
map. In the [0.5-2] keV energy band the X-ray surface bright- 
ness depends weakly on the gas temperature; for this reason the 
de-projected density profile using this band has a weak depen- 
dence on temperature variations. Inside rsoo the temperature can 
be fit with a constant (i.e. isothermal) value, at 10.4 +1.4 keV, 
marked by the horizontal dashed line in Fig|5] 

Since no "de-projected" temperature profile for Abell 2163 
from X-ray spectroscopic measurements is available in the lit- 
erature, we compute the mean weighted value of the gas tem- 
perature along the line of sight to compare our results with the 
published X-ray temperature values. The projected gas temper- 
ature is computed as Tproj = / WTdV / J WdV, where T is the 
de-projected gas temperature obtained from Abel inversion, and 
W is the weight function. As expected, the effect of projection 
is small compared to the eiTors in our temperature profile. We 
compare two different weighing schemes: the standard emission 
weight with W = n^MJ), and the weighing for a "spectroscopic- 
like" temperature as discussed by Mazzotta et al. (2004), us- 
ing W - rP'T^^^*. The difference between these two weighing 
schemes is negligible, as can be expected from the slowly vary- 
ing temperature profile of A2163 which shows no strong non- 
isothermal features. 

In Figure |6] we compare the emission-weighted temperature 
with published spectroscopic measurements from XMM-Newton 
and Chandra observations. The XMM-Newton temperature pro- 
file indicates a drop in temperature in the central region, which 
may be the remnant of the cold core of the original more massive 
cluster in this merging system. Gas temperatures derived from 
the APEX-SZ measurement are consistent with this feature. In 
the outskirts, gradually decreasing temperatures with incresing 
radii have been observed for relaxed clusters (e.g. Pratt et al. 
2007), and is expected from both theory (e.g. Frenk et al. 1999) 
and numerical simulations of cluster models (e.g. Roncarelli et 
al. 2006, Hallman et al. 2007). Although the large uncertainties 
in the present temperature profile beyond rsoo do not allow any 
conclusions as to a decrease in temperature at large radii, the 
data are fully consistent with such a trend. 



1000 1500 
r (kpc) 



Annular distance (arcminutes) 
4 6 8 




r (kpc) 

Fig. 5. De-projected radial density and temperature profiles of 
Abell 2163 from X-ray (XMM-Newton) and APEX-SZ data. 
Top: De-projected radial density profile (Icr confidence region). 
The dashed line is the X-ray derived spherical fi profile {j3 - 0.64 
and 0c - 91"), normalized to the derived central electron density 
value. Bottom: De-projected temperature profile with Icr uncer- 
tainties. The horizontal dashed line shows the best-fit isothermal 
value of 10.4 +1.4 keV within rsoo (~7.6')- Vertical arrows in 
both plots mark the estimated values of r2oo and r5oo- 



5.3.2. Gas mass and total mass 

The results of the mass analysis are shown in Figure [T] We have 
used median smoothing to the temperature and density profiles 
before computing their derivatives, and limited this analysis to 
within rsoo (approximately 1500 kpc) to avoid the large uncer- 
tainties at the outer radii. The errors in the total mass and gas 
mass fraction profiles reflect the combined random and system- 
atic errors in our analysis. We find that the cumulative mass pro- 
file obtained under the assumption of HSE with spherical sym- 
metry is in good agreement with the total mass obtained from 
weaklensing (Squires et al. 2007). The present analysis provides 
a better constraint on the total mass profile if the assumption of 
HSE is valid. The cumulative /gas values near the cluster cen- 
ter obtained from the non-parametric modeling are systemati- 
cally higher than the isothermal y6-model predictions. The drop 
to zero of the gas mass fraction at small radii is a well-known 
artefact of the isothermal yS-modeling of cluster ICM, and our 
non-parametric profile may be indicative of the actual ratio be- 
tween gas and total gravitational matter near the cluster center. 
At a radius of 1 Mpc and beyond, the /gas values from both meth- 
ods agree within Icr. 
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Fig. 6. Comparison of the projected radial temperature pro- 
file for Abell 2163 with spectroscopic X-ray measurements. 
The shaded region shows the emission-weighted temperature 
computed from the de-projected temperature and density pro- 
files. The squares and crosses are X-ray spectroscopic measure- 
ments from Chandra (Markevitch & Vikhlinin 2001) and XMM- 
Newton (Snowden et al. 2008), respectively. 



The resulting non-parametric gas mass fraction profile for 
Abell 2163 is close to the universal baryon fraction obtained 
from the WMAP 5-year data (Dunkley et al. 2009), and is also 
consistent with the X-ray analysis of Squires et al. (1997) for this 
cluster. The joint SZE/X-ray analysis of the gas fraction in clus- 
ters using OVRO/BIMA data (LaRoque et al. 2006) also yields 
a value of /gas in the range 0.15 - 0.16 for this cluster, using a 
non-isothermal double j6-model. The latter measurement is lim- 
ited to within r25oo, and at this radius the SZE-only isothermal 
value is approximately 30% lower than the double j6-model re- 
sult, consistent with the difference found in our analysis. 

It should be noted that the plotted mass and gas mass frac- 
tions are cumulative functions, and hence their values in the 
outer bins are correlated with the data in bins closer to the center. 
Furthermore, the errors in /gas are derived from the errors in the 
total mass and gas mass profiles, which are not independent. For 
these reasons, the errors in /gas are also correlated. 

6. Constraints from the SZE spectrum 

To compare the SZE decrement/increment values at different fre- 
quencies, we write them in terms of the relative change to the 
background CMB intensity. 



A/ = /o h{x) Qbeam 



Ar 

T^CMB 



(15) 



where x = hv/kTcuB, lo = ^ikT cmb)^ /(hc)^, f^beam is the beam 
equivalent solid angle in steradians, and h(x) relates the fre- 
quency dependence of the SZE expressed in temperature and in 
intensity as 



h(x) 



x^e^ 



(e' - 1)2' 

The total change of intensity is written as the sum of the thermal 
and kinematic components of the SZE as 



A/ = A/t 



A/k, 



where 



A/t = /(, y h(x) /(x, T,) 



(16) 
(17) 



10'=. 



10'". 



10' 



-r 






^ ^ 31- 
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Fig. 7. Top: Gas mass and total mass in Abell 2163 within rsoo 
(~1500 kpc), obtained from the de-projected density and tem- 
perature profiles. The shaded region shows the total mass within 
68% confidence levels. The dashed line is the gas mass, con- 
strained primarily from the X-ray data. The square boxes with 
error bars are the weak-lensing mass profile from Squires et 
al. (1997). Bottom: The gas mass fraction obtained from the 
isothermal analysis (solid line) and from the non-parametric de- 
projection (shaded region, 68% CL). The triangles with error 
bars are the results from the X-ray analysis of Squires et al. 
(1997), and the horizontal dot-dashed line represents the cos- 
mic baryon fraction from the WMAP 5-year result (Dunkley et 
al. 2009). 



and 



V,- 



(18) 



with x'r the radial (line-of-sight) peculiar velocity, which is posi- 
tive for a receding cluster. 

Relativistic corrections to third order in k-RTelnifC- are in- 
cluded in f{x, Te) according to Sazonov & Sunyaev (1998) (Note 
that the quantity h(x) fix, T^) is often denoted g{x, T^) in the 
literature). Higher order relativistic corrections are negligible 
given the precision of the current measurements. For the pur- 
pose of fitting the SZE spectrum, the Comptonization parameter 
y (defined in equation O is parameterized in terms of its central 
value Jo - This allows a simple model where the finite resolutions 
of different experiments are accounted for. All the spectral data 
are obtained with a radial fit using spherical /? models, and ac- 
counting for beam dilution. For APEX-SZ and LABOCA, we 
use the central decrement/increment values, corrected for dust 
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Table 3. SZE decrement/increment measurements used for the 
spectral fit. 



Wavelength (mm) 


Instrument 


A/ (MJy sr-') 


10 


OVRO/BIMA' 


-0.043±0.005 * 


2.1 


SuZIE^ 


-0.342±0.033 ^•'^ 


2.0 


APEX-SZ 


-0.317±0.035 


1.4 


SuZIE- 


-0.093±0.069 


1.1 


SuZIE- 


0.266±0.095 


0.86 


LABOCA 


0.633+0.094 



' LaRoque et al. (2002) 
^Holzapfeletal. (1997) 

3 



dust-corrected intensities from LaRoque et al. (2002) 
calibration corrected according to Hill et al. (2009) 




1000 2000 
BULK VELOCITY (km/s) 



contamination, derived in section |4] The data used for the fit to 
the SZE spectrum are given in Table [3] 

Due to the degeneracies between velocity, temperature and 
Comptonization, the present data is not sufficient for a simulta- 
neous constraint of all three parameters. Instead, we fix the ICM 
temperature to 10.4 keV from the joint X-ray and SZE analysis 
of section|5] and perform a least squares fit in yo and vv. We also 
fit these parameters using alternative values of the temperature 
in the range 8-14 keV. 

To estimate the errors in the fitted parameters, we perform a 
Monte Carlo simulation, in which 10,000 artificial data sets are 
created from the actual spectral data, adding random Gaussian 
offsets with the amplitude of the statistical noise at each fre- 
quency. To model the effect of the systematics, we also add 
random offsets from the most important systematic components 
from section l4~4l with the proper correlations between the differ- 
ent frequencies taken into account. Dust and primary CMB sig- 
nals are scaled to the estimated levels derived in section|4]at 150 
and 345 GHz and to the levels derived by LaRoque et al. (2002) 
at the other frequencies. To account for primary calibration un- 
certainties, we assume a 5% calibration uncertainty across the 
spectrum. 

The resulting error estimates, including systematics, are 
shown in Figure [8] The fitted radial velocity of the cluster is 
consistent with zero, in agreement with and marginally im- 
proving the constraints of LaRoque et al. Specifically, we find 
Vr = -140 ± 460 km/s and yo = 3.42 ± 0.32 x lO"'*, excluding 
systematic effects. The fit is dominated by random noise; we find 
that including the systematic errors in the modeling increases the 
error in yo by 12% and in by 18%. The SZE spectra resulting 
from the various fits are indicated in Figure |9l 

7. Conclusions 

1. We present SZE maps of the galaxy cluster Abell 2163 at 
two frequencies, showing the SZE decrement at 150 GHz 
from observations with the APEX-SZ bolometer camera, and 
the SZE increment at 345 GHz from observations with the 
LABOCA bolometer camera. The 345 GHz measurement is 
the highest resolution SZE image for this cluster to date, and 
the first large-area (> 10') imaging of a galaxy cluster at sub- 
mm wavelengths. 

2. An isothermal modeling of the SZE yields results consistent 
with the X-ray derived isothermal model fits, implying that 
the large-scale properties of the cluster, under the assumption 
of spherical symmetry, are well represented by such models. 

3. Using the APEX-SZ map in conjunction with XMM-Newton 
X-ray data we derive the (de-projected) radial density and 




1000 -500 500 

BULK VELOCITY (km/s) 



Fig. 8. Constraints in the Vr - yo parameter space of Abell 2163 
with priors on the ICM temperature, including systematic un- 
certainties. Top: solid contours indicate the 68% and 95% confi- 
dence regions using all available data and assuming kTf, - 10.4 
keV. The dashed contours correspond to excluding the present 
(APEX-SZ and LABOCA) measurements. Bottom: 68% confi- 
dence regions for a range of ICM temperatures (indicated) using 
all available data. 



temperature profile of the ICM under the assumption of 
spherical symmetry. The projected gas temperature profile 
is found to be fully consistent with previous X-ray spectro- 
scopic measurements of Markevitch & Vikhlinin (2001) us- 
ing Chandra data, and Snowden et al. (2008) using XMM- 
Newton data. 

The total mass profile of the cluster is obtained under the 
assumption of hydrostatic equilibrium. The resulting profile 
agrees well with the weak-Iensing mass profiles from previ- 
ous works (Squires et al. 1997, Radovich et al. 2008), and 
extends the profile out to r^oo- 

The gas mass fraction of the ICM obtained from our non- 
parametric analysis are consistent with the previous X-ray 
mesurements (Squires et al. 1997) and SZE/X-ray joint anal- 
ysis with double y6-modeI (LaRoque et al. 2006). We do not 
see any significant trend for increasing ICM baryonic frac- 
tion from r25oo to rsoo of the cluster. 

Using isothermal fits to the LABOCA and APEX-SZ mea- 
surements, we constrain the line-of-sight peculiar velocity of 
the cluster and the central optical depth to inverse Compton 
scattering, using the temperature obtained from the de- 
projection analysis. We find a peculiar velocity vv = -140 ± 
460 km/s, consistent with zero, and a central Comptonization 
yo = 3.42 + 0.32 X 10-4. 
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Fig. 9. SZE spectrum of Abell 2163 (points) and best- fit models 
using different priors on the ICM temperature: 8 keV (solid line), 
10 keV (long-dashed line), 12 keV (short-dashed line) and 14 
keV (dotted Une). 
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